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Abstract. We consider statistical analysis of double couple (DC) earthquake focal mech- 

^ ; 

^ anism orientation. The symmetry of DC changes with its geometrical properties, and the 

O 

^ number of 3-D rotations one DC source can be transformed into another depends on its 

symmetry. Four rotations exist in a general case of DC with the nodal-plane ambiguity, two 
J>~>" transformations if the fault plane is known, and one rotation if the sides of the fault plane are 

^ known. The symmetry of rotated objects is extensively analyzed in statistical material tex- 

ture studies, and we apply their results to analyzing DC orientation. We consider theoretical 
probability distributions which can be used to approximate observational patterns of focal 
mechanisms. Uniform random rotation distributions for various DC sources are discussed, 
as well as two non-uniform distributions: the rotational Cauchy and von Mises-Fisher. We 
discuss how parameters of these rotations can be estimated by a statistical analysis of earth- 
quake source properties in global seismicity. We also show how earthquake focal mechanism 
orientations can be displayed on the Rodrigues vector space. 
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1 Introduction 



This paper addresses two problems: the random rotation of double-couple (DC) earthquake 
sources and how symmetry properties of these sources influence their rotation angle distri- 
bution and their display. Properties of earthquake focal mechanisms and methods for their 
determination are considered by Snoke (2003) and Gasperini & Vannucci (2003). Ekstrom et 
al. (2005, and references therein) discuss their extensive work on evaluating seismic moment 
tensors for global earthquakes. 

In this paper we considered only the double-couple earthquake focal mechanism. For 
tectonic events non-double-couple mechanisms like the CLVD are likely due to various sys- 
tematic and random errors in determining the mechanism (Prohlich & Davis 1999; Kagan 
2003, 2009). These results suggest that routinely determined CLVD values would not reliably 
show the deviation of earthquake focal mechanisms from a standard DC model. 

Snoke (2003) and Gasperini & Vannucci (2003) consider several equivalent representations 
for double-couple sources and their properties, and provide mathematical expressions for 
their mutual transformation. Krieger and Heimann (2012, and references therein) review 
routines for plotting moment tensors and focal mechanisms. 

Two general techniques can be employed to study the 3-D rotation: orthonormal rotation 
matrices and normalized (unit) quaternions. The quaternion method has been used to 
evaluate these rotations in many investigations of earthquake focal mechanisms (see, for 
example, Kagan 1991; Prohlich & Davis 1999; Kagan 2009; Kagan & Jackson 2011). Kagan 
(2007) explains how 'ordinary' matrices and vectors can be used to obtain 3-D rotation 
parameters. Below we comment on the advantages and drawbacks of both methods. 

Altmann's (1986) book was a first monograph specifically dedicated to 3-D rotations 
[group 5*0(3)] and quaternions. At present quaternions are widely used to describe rota- 
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tions in space satellite and airplane dynamics (Kuipers 1999) and simulations of virtual 
reality, robotics and automation (Hanson 2006; Dunn & Parberry 2011). These last three 
monographs explain quaternions in a more accessible manner. Many journal articles (see ref- 
erences in these monographs and in Kagan 2009) discuss practical application of quaternions 
for analyzing the 3-D rotations. 

However, the above publications do not consider the symmetry properties of rotated ob- 
jects or how symmetry influences orientation analysis. As Kagan (1990; 1991; 2007; 2009) 
indicated, the techniques considered in those publications cannot be used for the DC source 
orientation studies without major modifications, because of the DC symmetry properties, 
described in Section [2l The only scientific discipline where symmetry is extensively consid- 
ered in 3-D rotation analysis is study of material texture (Handscomb 1958; Mackenzie 1958; 
Grimmer 1979; Frank 1988; Heinz & Neumann 1991; Morawiec, 2004; Meister & Schaeben 
2005; Schaeben 2010). In this paper (Section [3]) we discuss applying these results to DC 
source investigations. In Section H] we consider theoretical probability distributions used 
to approximate observational patterns of focal mechanisms. Section O is dedicated to the 
statistical analysis of earthquake source properties in global seismicity. Using the results of 
Section [3] we also show how earthquake focal mechanism orientations can be displayed in the 
Rodrigues vector space. Section E] and Section [7] summarize our results. 



2 Focal mechanism symmetry 

Depending on the known properties of double-couple earthquake focal mechanism, we con- 
sider three types of earthquake source symmetry (Kagan 1990): 

• 1. DCl - double couple with no symmetry or the identity (J) symmetry, if the focal plane 
and its sides are known; 
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• 2. DC2 - double couple with C2, order 2 cyclic symmetry, i.e., the focal plane is known, 
but its sides are not; 

• 3. DCA - double couple with nodal-planes that not distinguishable; it has D2, order 2 
dihedral symmetry. 

These earthquake source symmetries correspond to the following crystallographic sym- 
metries considered in material texture analysis (see, for instance, Morawiec 2004). DC4 has 
an orthorhombic symmetry (as in a rectangular right parallelepiped or a rectangular box 
with unequal sides); DC2 has a monocline symmetry (as in a 3-D prism with two angles of 
90° and one arbitrary angle); DCl has a triclinic, or no symmetry. 

Fig. [1] displays the geometry of the DC source (Aki & Richards 2002). It represents the 
quadrupolar 'beachball' radiation patterns of earthquakes. The focal plots involve painting 
on a sphere the sense of the first motion of the primary P-waves: solid for compressional 
motion and open for dilatational. The two orthogonal nodal planes separating these areas 
are the fault and the auxiliary planes. During routine determination of focal mechanisms, it 
is impossible to distinguish between these planes, a property called 'nodal-plane ambiguity.' 
The planes' intersection is the null-axis (called b-axis), the p-axis is in the middle of the 
open lune, and the t-axis is in the middle of the closed lune. These three axes are called 
the 'principal axes of an earthquake focal mechanism,' and their orientation defines the 
mechanism. 

To make the focal mechanism picture unique, the eigenvectors are pointed down in regular 
representations. However, the handedness of the coordinate system formed by the vectors 
can change as the result of such an assignment. The systems of the opposing handedness 
cannot be rotated one into another. In most of our considerations, we use the right-handed 
coordinate system placed at each earthquake centroid. 

Fig. |2] displays four examples of the right-handed coordinate system for a DCA source. 
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The system can be arbitrarily rotated, and the handedness of the system is preserved. The 
left-handed system can be obtained in this picture if one inverts the direction of any indi- 
vidual axis or of all three axes. If the direction of two axes is reversed, the handedness of 
the system is preserved. 

The earthquake focal plane can often be determined by inverting the higher-rank point 
seismic moment tensors (McGuire et al. 2001; Chen et al. 2010) or by the aftershock pattern. 
The face/side (up/down or foot /hanging wall) of a focal plane generally is unknown. In such 
a case the DC focal mechanism has a C2 symmetry; we call it DC2. 

Finally, the face or the side of the focal plane can be known; it is shown in Fig. [T] by 
symbols A and B. Such a source is called DCl. Fig. |3] illustrates the difference between 
DC2 and DCl mechanisms. If a cylinder of one material is rotated in a half-space of another 
material, the two DCl sources '1' and '2' would have a rotation angle of 180°. If they are 
considered as DC2 sources, two angles 0° and 180° are possible to rotate one mechanism 
into another. 

In Fig. [21 any of the configurations appropriately rotated can represent a DCl source, but 
only pairs (a)-(b) or (c)-(d) correspond to a DC2 mechanism. All four diagrams correspond 
to a DCA source. 

3 Earthquake focal mechanism and material texture 
statistics 

Frank (1988) proposes using the Rodrigues vector space to represent 3-D rotation of sym- 
metrical objects. This representation has an advantage: under any transformation of the 
Rodrigues map corresponding to a change of the reference orientation, straight lines trans- 
form into straight lines, and planes into planes. 
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For an object with the non-identity symmetry, accepted points he in a region around 
the origin, which is caUed 'the fundamental zone of the map.' It is a polyhedron, bounded 
by the planes which are orientationally equidistant between the origin and the neighboring 
equivalent point by a symmetry rotation to the origin. Any points lying outside one of these 
planes have an equivalent point lying inside the fundamental zone (Frank 1988). For an 
orthorhombic crystal with three orthogonal axes, the fundamental zone is a cube, with its 
six faces orthogonal to the axes at a distance from the origin of tan7r/4 = 1 (Frank 1988). 
The cube is surrounded by three neighboring zones, each divided into two at infinity. 

Fig. m shows the fundamental zone for a DCA source (Heinz & Neumann, 1991, Fig. 7). 
It is a cube with corner coordinates xi = ±1.0; X2 = ±1.0; X3 = ±1.0. Owing to the DCA 
symmetry, an octant of the cube contains full information about the orientation distribution 
for uniformly random rotation. This octant is called the 'MacKenzie cell' (Morawiec & Field 
1996, see also their Fig. 1 displaying the cells for the and -D4 symmetries). 

Each point inside the cube uniquely corresponds to a certain orientation/rotation with 
a minimum rotation angle $min ^ 120°. The points in the inscribed sphere of the cube 
correspond to the rotations with angles $niin < 90°. The other three rotation angles 
are situated outside the fundamental zone. For example, the point of the zero rota- 
tion is located at the cube center, whereas three other rotation points are at infinity: 
Xi = ±00; X2 = ±00; X3 = ±00, corresponding to 180° rotations. These points at ±00 
are equivalent (Altmann 1986). Similarly, for any point inside the cube, three points outside 
correspond to the rotations with angles $ > $min- 

However, when a point moving orthogonally from the origin reaches a cube face, it 
simultaneously appears on the opposite face: two 90° rotations produce the same effect 
(Fig. in]). This means that when we determine the minimum angle $min for cyan point 
rotations shown in Fig. [5] using the program developed by Kagan (1991), we find two equal 
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solutions. The remaining two angles are greater than $min (see also the next Section). 

If a point on one face moves to an edge, the 'identical' point on the opposite face si- 
multaneously moves to another edge until both points reach the middle of the edges. This 
orientation corresponds to the rotation $ ^ 109.5°, and Fig. [6] shows that there are three 
equivalent points at the edges. The third point appears as it moves from the outside of 
the cube to the third edge. As in Fig. [5], this means that three equal angles $min would 
be obtained. Finally, when a point is at a vertex, as shown in Fig. [TJ three other vertices 
correspond to the same rotation $ = 120° (Frank 1988), i.e., all four rotation angles are 

^min- 

This arrangement of the orientations for the rotation angles $ > 90° describes a complex 
topology for DC4 source rotation. This topology involves projective or Mobius transforma- 
tion (Altmann 1986; Frank 1988). Full analysis of the DC4 source orientation, when and if 
performed, would involve very intricate investigations of rotation angle transformations due 
to source symmetry. 

The Rodrigues space has no special advantages in displaying the orientation distribution 
for two other sources: DC2 and DCl. For DC2 the fundamental zone overlaps the entire 
Rodrigues space. For C2 symmetry the fundamental zone is bounded by two planes perpen- 
dicular to the b-axis, each at the distance tan(7r/4) = 1 from 0. For a DCl source, the whole 
Rodrigues space up to infinity is included. In these cases other spaces are more convenient 
in displaying the 3-D rotation distribution (Frank 1988; Morawiec & Field, 1996). Altmann 
(1986, pp. 164-176) explains the projective or Mobius topology of rotations in the quaternion 
parametric ball for non-symmetrical objects. 
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4 Rotation of double couple earthquake sources 



4.1 DC symmetry and rotation angle 

In our earlier paper (Kagan 1991) we considered the inverse problem of the DC4 source 
rotation, i.e., given two earthquake focal mechanisms to determine all the 3-D rotations by 
which one mechanism can be rotated into another. As we showed, the DCA symmetry results 
in four such rotations with angle $ range 0° < $ < 180° (cf. Fig. For most practical 
purposes, the rotation with a minimum angle $min can be selected. Following the material 
texture designation (Morawiec 2004, p. 115), we sometimes refer to general rotation angles 
as 'misorientation angles' and the minimum rotation angle as the 'disorientation angle.' 

Here we consider how to compute the rotation of a DCl source needed to align it with 
some reference DCl source. It is unlikely that sufficient data would exist on DCl sources 
for a statistical study of their distribution. However, in some cases we need to measure their 
angle of disorientation or the angular distance. Moreover, if the disorientation angle between 
two DCl sources is large, it would be almost impossible to identify their fault planes and 
plane faces. Small rotation angles for such focal mechanisms are likely to correspond to $min 
for the DCA source, as mentioned above. 

To compute the disorientation of a DCl source, we can modify our FORTRAN program 
listed in Kagan (1991). When a fault plane and its face are known, a focal mechanism would 
be better specified through a fault plane geometry (Aki & Richards 2002, Figs. 4.13 and 
4.20) with three angles: strike or azimuth (0), dip (S), and rake (A). Usually the range of 
these angles is taken as follows: 0° < < 360°, 0° < 5 < 90°, -180° < A < 180°. The 
problem arises when comparing two sources if the dip (6) of one focal plane exceeds 90°, 
so that a foot wall of one mechanism becomes a hanging wall for another source (Aki & 
Richards 2002). To simplify the calculations in our program (see below), we extend the 6 
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range to 180°. 

If the face/side of a fault plane is unknown, as shown in Fig. [3l we need to calculate 
the second angle of the rotation for the DC2 source. As with the DCl source, the data on 
DC2 sources are sparse and insufficient for a statistical study, but we need a technique to 
measure their angles of disorientation. An easy way to accomplish this measurement would 
be to change the strike of the fault plane by 180° and change the rake sign. The modified 



DClROT.FOR program is available at http://jumpy.igpp.ucla.edu/~kagan/dclrot.for 



In this program we use the quaternion technique to determine the rotation angle and the 
rotation axis parameters to transform one DC source into another. Quaternions are used 
because for rotation angle $ close or equal to 180°, the matrix method cannot determine the 
rotation axis parameters (Kagan, 2007). 

Kagan (2009, Appendices A and B) discusses normalized quaternions and their relation 
to the DCA source. Representations for the DCA by seismic moment tensors as well as by 
orthonormal matrices are considered. A quaternion representation allows a relatively easy 
determination of the rotation angle and the axis parameters for the DCA earthquake focal 
mechanism (Kagan 1991). 

If only the rotation angle is needed, then one can use a scalar (dot) product of two 
quaternions (Hanson 2006, p. 65; Dunn & Parberry 2011, p. 255) to determine the angle: 

cos($/2) = q'^ ■ q" = q\ + q^, + + ql q\ , (1) 

where q* are normalized quaternions for each DC source and g-,- are the quaternion's com- 
ponents. 

In our program we first compute the orthonormal matrix for each DC source and then 
determine the corresponding normalized quaternion. There is a possibility of losing precision 
when converting a matrix to a quaternion (Shepperd 1978; Horn 1987). A certain computa- 
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tion technique should be apphed to avoid this. We used a similar technique in our programs 
DCROT.FOR (Kagan 1991, the end of the SUBROUTINE QUATFPS) and in dcIrot.for (see 
above) . 

4.2 Rotation angle distribution 

4.2.1 Uniform random rotation of DC sources 

The distribution of the uniform random rotation for DC sources constitutes a reference for 
distributions occurring in earthquake focal mechanisms where we expect the distributions to 
be partially random. These stochastic distributions can be analytically calculated by taking 
into account the sources' symmetry. 

A distribution of the minimum angle $min for a uniform random rotation of the DC4 
source was obtained by Kagan (1990, Eqs. 3.1-3.3), using the results by Handscomb (1958) 
and Mackenzie (1958) for the random disorientation of two cubes. The probability density 
function (PDF) is 



/($) 



(4/7r)(l - cos$) for 0<$<7r/2; 



(2) 



/(*) 



(4/7r)(3sin$ + 2cos$ - 2) for 7r/2 < $ < $5 ; 



(3) 



and 





(G/tt) 2 sin$ arccos 



f l+cos <!> \ 
V -2 cos * ) 



1/2 



(1 — cos$) arccos 



l+cos ^ 
—2 cos $ 




for $5 < $ < 1^ , 



(4) 



where 




— arccos 
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For the DC2 source a similar PDF is 



/($) = (2/7r)[l - cos($)] for < $ < 7r/2 ; 



(6) 



and 



/($) = (2/7r) sin($) for 7r/2 < $ < tt . 



(7) 



For the DCl source the function is 



/($) = (I/tt) [1 - cos($)] for < $ < TT . 



(8) 



Grimmer (1979), also following Handscomb (1958) and Mackenzie (1958) results, ob- 
tained similar analytic expressions for a completely random rotation of orthorhombic, mon- 
oclinic, and triclinic crystals (equivalent in symmetry to the DC earthquake source with 
various restrictions described above). He listed median angles as well as their mean and 
standard deviations for all these distributions. Morawiec (1995; 2004, pp. 117-119) derived 
these distributions by integration in the Rodrigues space. 

4.2.2 Non-uniform distributions of random rotations 

Two non-uniform rotation angle distributions are useful in analyzing earthquake focal mecha- 
nism rotation: the rotational Cauchy law (Kagan 1982, 1992) and von Mises-Fisher/Bingham 
rotational distribution (Kagan 1992, 2000; Schaeben 1996; Mardia & Jupp 2000, pp. 289-292; 
Morawiec 2004, pp. 88-89). 

The Cauchy distribution is especially important for representing earthquake geometry, 
since it can be shown by theoretical arguments (Zolotarev 1986, pp. 45-46; Kagan, 1990) and 
simulations (Kagan 1990) that the stress tensor in a medium with random defects follows 
this distribution. The Cauchy law is a stable distribution (Zolotarev, 1986). The stable 
distributions are essential for two reasons: a) They are invariant under addition of random 
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variables; b) Stable distributions have a power-law tail, i.e., they are asymptotically scale- 
invariant. 

The probability density function (PDF) of the rotational Cauchy distribution can be 
written as (Kagan 1982; 1990) 



where A = tan($/2). The scale parameter k of the Cauchy distribution represents the 
degree of incoherence or complexity in a set of earthquake focal mechanisms. The cumulative 
rotational Cauchy distribution can be written as 



The Cauchy distribution is assumed to be axisymmetric on the quaternion hypersphere 
S^. This means that the rotation axis poles are distributed uniformly over a regular S"^ 
sphere. For a general case, the axes distribution for earthquake focal mechanisms may need 
to be specified as non-uniform. In that case certain rotations would be preferred depending 
on the focal mechanism of a reference event. However, we have not yet advanced to this 
stage (see Section [5]). 

The von Mises-Fisher/Bingham distribution for the 3-D orientation is widely discussed 
in literature (Schaeben 1996; Mardia & Jupp 2000; Morawiec 2004). Schaeben (1996) and 
Morawiec (2004) show that this distribution is essentially equivalent to the Bingham dis- 
tribution. The von Mises-Fisher distribution is a Gaussian-shaped function concentrated 
near the zero angle. This distribution can be implemented to model random errors in de- 
termining focal mechanisms. The distribution has many forms. However, even the simplest 
axially symmetric expressions, due to the complexity of normalization, represent difficult 
computation. 




for 180° >$ > 



(9) 




(10) 
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For small values of the standard error a^, the von Mises-Fisher-type distribution is equiv- 
alent to the rotational Maxwell law used by Kagan & Knopoff (1985) and Kagan (1992). 
The latter distribution is obtained by generating a 3-D normally distributed random variable 
u {ui, U2, U3) with the standard deviation (cr^^, cr^j, c^g) and then calculating the unit 
quaternion 

(Jo = l/\Jl + ul + ul + ul, 
Qi = Ui/^l + ul + u'^ + ul, fori = 1,2, 3. (11) 

The 3-D rotation angle is calculated 

$ = 2arccos(go) ~ 2 Q3:ccos(i — (uj + ul + ul)/2^ 



2 arcsin 2^Jul + U2 + ul. (12) 

The final expression is twice the length of a vector in the 3-D space. 

Since components of vector u are normally distributed, the angle $ (in degrees) follows 
the Maxwell distribution with 

(7$ = 360(7„/7r, (13) 

where we assume that all components of Uu are equal (uu^ = — cr^^ — Uu)- For 180° > 
$ > 0° the Maxwell PDF is 

/2" $2 

^($) = J-x^xexp -$V(2cr|) ■ (14) 

V TT (7$ L J 

This equation describes the distribution of a vector length in three dimensions, if the vector 
components have a Gaussian distribution with a zero mean and a standard error The 
Maxwell cumulative distribution function (CDF) is 

v]>($) = erf - J- X — X exp [-$7(24)1 , (15) 
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where erf (.) is an error function. 

The major problem with these non-uniform random distributions is that they do not 
consider the symmetry of the rotated object. When rotation angles are small, the distribution 
is concentrated around the zero angle neighborhood. For the DCA source as shown in Fig. HI 
almost all distribution density would be inside the fundamental zone. However, for more 
spread out angle distributions, we should account for cases where the rotation angle exceeds 
the maximum angles (see, for example, Eqs. |2]|5]). Then the distribution would be folded 
back into the fundamental zone. 

Mason & Schuh (2009) propose to convolve angle distribution with appropriate 3-D 
spherical or 4-D hyperspherical harmonics to obtain a new angle distribution which fits into 
the fundamental zone. It is not clear whether such calculations can be made analytically. 
Simulation seems the only practical way to transform both Cauchy and von Mises/Fisher 
distributions for the D2 symmetric case (i.e., for maximum $min rotation angle 120°). Kagan 
(1992, Fig. 3c) produced such distributions. Fig. |8] below plots the appropriate Cauchy 
distribution which is reduced to ^rnm < 120°. 

5 Focal mechanisms statistics 
5.1 Disorientation angle statistics 

Since there is no general model of earthquake focal mechanism distribution, we need to study 
the distribution of mechanisms in earthquake catalogs empirically to infer their properties. 
How various tectonic and geometrical factors shape the distribution of earthquake sources 
should be studied as well. Such investigations are difficult because we are dealing with a 
multidimensional stochastic point process: earthquake size, occurrence time, location, and 
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source parameters serve as potential inputs to the distributions. 

In this paper we are mostly interested in the distributions of rotation between two earth- 
quake focal mechanisms. Even if we fix earthquake time, space, and magnitude interval, the 
DC rotation distribution depends on at least three variables: the rotation angle and two 
spherical coordinates of a rotation axis pole. Displaying all three degrees of freedom in a 
distribution presents a difficult problem. Therefore, in our previous investigations we stud- 
ied partial distributions. For example, Kagan (1992, Figs. 6-9; 2009, Figs. 9,10) obtained 
various distributions of the rotation angle $min between two focal mechanisms. Below we 
first update our most important results on the distribution of the rotation angle $min, and 
then analyze a three-dimensional distribution of rotation angle and the axes in the Rodrigues 
space. 

We used the Global Centroid Moment Tensor catalog (Ekstrom et al. 2005), referred 
to subsequently as GCMT. This catalog employs relatively consistent methods and reports 
tensor focal mechanisms. The GCMT catalog started in 1977, and is complete only for 
earthquakes with magnitudes of about 5.8 and larger. The present catalog contains more 
than 36,000 earthquake entries from 1977/1/1 to 2011/12/31. 

Fig. |8] displays cumulative distributions of the rotation angle $min for shallow earthquake 
pairs with the magnitude threshold = 5.0 that are separated by a distance of less than 
50 km. We study whether the rotation of focal mechanisms depends on where the second 
earthquake of a pair is situated with regard to the first event. Thus, we measure the rotation 
angle for centroids located in 30° cones around each principal axis of the first event (see 
curves, marked the t-, p-, and b-axes). 

The curves in Fig. [8] are narrowly clustered, with about 95% of angles less than 90°, 
within an inscribed sphere of the fundamental zone (Fig. S]). This pattern can be compared 
to the uniform rotation (Eqs. |2]|4]) for which 72.7% (2 — 4/7r) of angles are within 90°. The 
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curves are obviously well approximated by the DC rotational Cauchy distribution (Eq. |9]). 
This distribution is characterized by a parameter a smaller K-value corresponds to the 
rotation angle $min concentrated closer to zero. Thus, regardless of spatial orientation, all 
earthquakes have focal mechanisms similar to a nearby event. Earthquakes in the cone 
around the b-axis correspond to a smaller K-value than events near the other axes. These 
results are similar to those shown in Fig. 6 by Kagan (1992) or Fig. 9 by Kagan (2009). 

Fig. [9] and Fig. [10] show the disorientation angle distribution for the magnitude cutoff 
rrit = 5.8. As may be expected for the higher magnitude, the angles are concentrated closer to 
zero, and the difference between the curves corresponding to various cones increases. Maxwell 
distribution curves are shown to illustrate possible behavior of the angle distributions near 
zero. The cr$ = 7.5° parameter of the distribution is small, compared to the distribution 
range (120° > $ > 0°). Therefore, the curves are concentrated close to zero; we do not need 
simulation to consider the curve behavior for large values of $min, as done, for example, in 
Fig. 3c by Kagan (1992). 

The difference in the distribution curves corresponding to various focal mechanism axes 
suggests that the Cauchy distribution parameter k depends upon the geometry of a fault 
system. Contrary to our assumptions (Eq. |9]), poles of rotation axes are not uniformly 
distributed over the sphere. 

5.2 Distributions of rotation axes 

Mackenzie (1964) derived the distribution of the rotation axes for cubic symmetry. Moraw- 
iec (1996) obtained distributions of rotation axes for any symmetric object encountered in 
material texture analysis. Using his results we can write down the distribution of rotation 
axes for the D2 symmetry: the DCA source. The Mackenzie cell is shown in Fig. |H We 
designate the coordinate 1, 2, 3, and the distribution depends on distance 
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from the origin. As seen in Fig. HJ the distribution should have a 3-fold cyclic symmetry C3 
around the origin or around the cube vertex. Then the PDF for the axes density is 



Fig. [m displays the p{p) density. Its values at p = 1, p = -\/2, and p = \/3 correspond 
to appropriate values for the D2 symmetry (Morawiec 1996, Table 2). These p- values cor- 
respond to the rotation angles $ = 90°, 109.5°, and 120.0°, respectively. For $ < 90° the 
rotation axes are distributed uniformly over the S'^ sphere, but they intersect the sphere 
near the cube vertex close to $ = 120°. 

In Fig. [12] we show the distribution of the rotation poles for the second earthquake focal 
mechanism on a reference sphere of the first event. Because of the symmetry of the DC4 
source, we reflect the point pattern on our reference sphere at the planes perpendicular to 
all axes. Thus, the distribution can be shown on a spherical octant. We use the Lambert 
azimuthal equal-area projection. The points concentrate near the projections of the far 
edges of the MacKenzie cell (see Fig. Hj) and around the cube vertex which corresponds to 
the disorientation angle $ = 120°. 

5.3 Rodrigues space statistics and display 

A major problem in the orientation visualizing is the high dimensionality of the 3-D rotation 
space: the orthogonal matrices are characterized by nine values, the seismic moment tensor 
requires five or six variables, and the normalized quaternion needs four values. The real 
number of degrees of freedom for a 3-D rotation is three. Thus, in principle, an orientation 
distribution can be shown in a 3-D diagram. 

Frank (1988), Neumann (1992) and Morawiec & Field (1996) propose using the Rodrigues 



P(p) 



- arctan(p) - p/(l p^) 



(16) 



where the distance p 




|. For small value of p, p (p) cx p^. 
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vector space to display the disorientation of symmetric objects in a fundamental zone as a 
point in the space. The point coordinates are calculated as follows: the length of a vector ( 
is 

C = tan($/2). (17) 
where $min is the rotation angle. Three coordinates of a point in the zone are 

Xi = ( X sin(6') sin(0); 
X2 = C X sin(6') cos(0); 

X3 = C X cos(^) , (18) 

where 6 is the colatitude, and is the azimuth of the rotation axis. As shown in Fig. HJ we 
identify Xi with the p-axis; similarly X2 = t and X3 = b. 

We obtain a distribution diagram for a set of disorientations. One way to display such 
a diagram of the fundamental zone is through stereo-pairs (Neumann 1992). Morawiec & 
Field (1996) display the distribution of misorientation parameters by points in some parallel 
sections of the fundamental zone. 

Fig. dSl shows a distribution for randomly rotated DCi sources in a central section of 
the fundamental zone of the Rodrigues space. Fig. [H] displays a similar distribution of the 
earthquake focal mechanism orientation in the GCMT catalog. As Figs. IHlfTOl demonstrate, 
the distribution of the rotation angles for earthquake sources is strongly concentrated close 
to 0°. If we compare Fig. [13] and Fig. [141 this concentration is marked in the fundamental 
zone display. 

Table [T] summarizes earthquake focal mechanism disorientation patterns in the funda- 
mental zone of DCA. The total number of events with the magnitude above the threshold 
is shown, as well as the total number of pairs Np with centroids separated by less than 50 km. 
Nc is the pair numbers in the central zone shown in Fig. [TH Although the central section 
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occupies only 5% of the zone, close to 50% of the pairs are there due to a tight concentration 
of rotation angles near the zero value. For simulated events in the central zone, the number 
Nc is about 7% of the total (see Fig. [T^ . 

We also display the correlation coefficients of the point scatter field. Whereas for earth- 
quake focal mechanisms the coefficients pbp and pu are close to zero, the ppt and p'p^ coefficients 
are non-zero, testifying to a certain pattern of focal mechanism rotation. Fig. 11 in Kagan 
(2009) also shows that rotation axes are concentrated closer to the t-axis. All correlation co- 
efficients for the simulated mechanisms are around zero. A more detailed statistical analysis 
of this pattern will be carried out in our future work. 

The values of the average rotation angle and its standard deviation ( $ ± cr$) show 
that for larger earthquakes both variables are smaller. This may be caused by a higher 
accuracy in determining focal mechanism for stronger shocks (Kagan 2003). For simulated 
focal mechanisms, the $ ± 0"$ values are close to the theoretical estimates for orthorhombic 
symmetry (Grimmer 1979). 

6 Discussion 

Quantitative study of earthquake focal mechanisms is an important prerequisite for under- 
standing earthquake rupture. Though these investigations began in the mid 1950s, pub- 
lications have been mostly descriptive until now; relatively little modelling and rigorous 
statistical analysis have been performed. A major difficulty in analyzing focal mechanisms 
is both the high dimensionality and non-commutativity of the 3-D rotations. This presents 
a major challenge in analyzing a set of focal mechanisms. 

Several papers (Kagan, 1992, 2000, 2009) have investigated statistical distributions of 
earthquake focal mechanisms. We found that the disorientation angle is close to zero for 
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spatially close earthquakes, and the angle decreases if the inter-earthquake time interval ap- 
proaches zero. We also showed (Kagan, 2009) that the CLVD component of focal mechanism 
tensor is either zero or close to zero for most geometric barriers proposed as common features 
in an earthquake fault system. 

However, the major challenges in describing and understanding the distributions of focal 
mechanisms still remain. As we see from Figs. IHlfTOl the angle distribution is not axially 
symmetric: in certain directions $mm is larger than in others. Thus, the distributions used 
to approximate the angle pattern, like the rotational Cauchy distribution (Eqs. [Ql fTOj) . need 
to be made more complex. 

The distribution of rotation axes was not investigated as thoroughly as that for disori- 
entation angles. There is still no theoretical model for approximating empirical data, but 
applying the Rodrigues space may render such analysis more manageable. 

However, even these limited results contribute significantly to understanding of earth- 
quake focal mechanism properties and allows certain quantitative applications for seismic 
risk evaluation. Kagan & Jackson (2011) explain that the forecasted tensor focal mechanism 
enables calculating an ensemble of seismograms for each point of interest on the Earth's 
surface. Moreover, the focal mechanism distribution allows us to estimate fault plane orien- 
tation for past earthquakes, through which we can identify a preferred rupture direction for 
future events. 

The angle $min has also been used to directly compare moment tensors from two different 
earthquakes (Okal et al. 2011). It has as well been applied in comparing moment tensors 
computed for the same events through different techniques (Frohlich & Davis 1999; Pondrelli 
et al. 2007; Yang et al. 2012). Such comparisons can help refine the moment tensor algorithms 
and lower their computational cost, since they can reveal the relative importance of various 
assumptions implicit in the algorithms. 
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7 Conclusions 



• 1. The symmetry properties of an earthquake double-couple focal mechanism are 
considered. Given available seismological and geologic information, three symmetries are 
possible: D2 (dihedral symmetry), C2 (cyclic symmetry), and / (identity). Determining the 
orientation or disorientation of the source depends on its symmetry. 

• 2. Quaternion representation is the most convenient tool for analysing a double-couple 
3-D rotation. A 3-D rotation requires at least three degrees of freedom for its characteriza- 
tion. 

• 3. Several theoretical distributions to describe a 3-D rotation of double-couples are 
presented: random rotation, rotational Cauchy, and von Mises-Fisher. 

• 4. The Rodrigues space, so extensively used in material texture analysis is applied 
to display and analyze earthquake focal mechanism distribution. This space allows us to 
represent 3-D patterns of how symmetric objects are oriented. 

• 5. We illustrate the proposed methods by statistically analyzing the GCMT catalog of 
earthquake focal mechanisms. 
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Tabic 1: Properties of disorientation point scatter in the fundamental zone of a DCi source 



# 


rrit 


N 






Ppt 


Pbp 


Pbt 


Ppt 




1 


5.6 


9,615 


43,611 


18,381 


0.126 


0.026 


0.013 


0.170 


31.6 ± 27.5 


2 


5.8 


6,160 


19,367 


8,725 


0.117 


0.022 


-0.001 


0.127 


29.9 ± 26.9 


3 


6.25 


2,154 


2,741 


1,244 


0.274 


-0.016 


-0.026 


0.257 


28.5 ± 26.1 


4 


Simul. 


25,000 


25,000 


1,754 


-0.012 


-0.002 


0.002 


-0.015 


75.0 ± 21.0 



Notes: The GCMT catalog time interval is 1977/1/1-2011/12/31. N is the total number 
of events with magnitude m > rrit; Np is the total number of event pairs; is the number 
of event pairs in the central section; pp^, pf,p, and pbt are the correlation coefficients for all 
points; p^^ is the correlation coefficient for all points within the central section; $ ± 0"$ is 
the average disorientation angle and its standard deviation. 
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Figure 1: 

Schematic (beachball) diagram of the DC earthquake focal mechanism and its quadrupole 
radiation patterns. The null (b) axis is orthogonal to the t- and p-axes, or it is located on 
the intersection of fault and auxiliary planes, i.e., perpendicular to the paper sheet in this 
display. The n-axis is normal to the fault-plane; u is a slip vector. 
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Figure 2: 

Four schematic diagrams of earthquake focal mechanism, having the DCA symmetry. The 
right-hand coordinate system is used. 
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Figure 3: Schematic diagram of earthquake focal mechanism DC2. 
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Figure 4: 

Fundamental zone display for DC A source. The b, p, and t-axes of the source are shown. 
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(a) 



(b) 





Figure 5: 

Fundamental zone display for DC 4 source. Colors show two face points corresponding to 
one source orientation with the angle $ > 90°. 
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Figure 6: 

Fundamental zone display for DCA source. Colors show four sets of three edge points 
corresponding to one source disorientation with the angle $ ~ 109.5°. Compare to Fig. [5] 
where two edge points are result of points moving on opposing faces. The third point appears 
as it moves from outside the fundamental cube to the third edge. Three other sets of points 
are similarly produced. 
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Figure 7: 

Fundamental zone display for DCi source. Colors show two sets of four vertex points 
corresponding to one source orientation with the angle $ = 120°. 
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CMT 1977-2011, shal-R-0— 50km, M5.0: o all, xT, + P, * B 




20 40 60 80 100 120 

Rotation Angle Degrees 

Figure 8: 

Cumulative distributions of rotation angles for pairs of focal mechanisms of shallow earth- 
quakes (depth 0-70 km) in the GCMT catalog 1977/01/01-2011/12/31; centroids are sepa- 
rated by distances between 0-50 km, magnitude threshold rrit = 5.0. The total number of 
events is 26,986. Lines from left to right: filled circles are centroids in 30° cones around the 
b-axis; dashed line is for the Cauchy rotation with k = 0.1; circles - all centroids; crosses - 
centroids in 30° cones around the p-axis; x-signs - centroids in 30° cones around the t-axis; 
right solid line is for the random rotation. 
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CMT 1 977-201 1 , shal-R-0— 50km,M5.8: o all, x T, + P, * B; — Cauchy; - Maxw. 
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Figure 9: 



Cumulative distributions of rotation angles for pairs of focal mechanisms of shallow earth- 
quakes (depth 0-70 km) in the GCMT catalog 1977/01/01-2011/12/31; centroids are sepa- 
rated by distances between 0-50 km, magnitude threshold mt = 5.8. The total number of 
events is 6,160. Lines from left to right: filled circles are centroids in 30° cones around the 
b-axis; circles - all centroids; crosses - centroids in 30° cones around the p-axis; x-signs - 
centroids in 30° cones around the t-axis; dashed line is for the Cauchy rotation CDF with 
K = 0.075; left solid line is for the Maxwell rotation CDF f|T5|) with cr$ = 7.5°; right solid 
line is for the random rotation. 
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CMT 1 977-201 1 , shal-R-0— 50km,M5.8: o all, x T, + P, * B; — Cauchy; - Maxw. 
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Figure 10: 



Histograms of rotation angles for pairs of focal mechanisms of shallow earthquakes (depth 0- 
70 km) in the GCMT catalog 1977/01/01-2011/12/31; centroids are separated by distances 
between 0-50 km, magnitude threshold rrit = 5.8. The total number of events is 6,160. Lines 
from left to right: filled circles are centroids in 30° cones around the b-axis; circles - all 
centroids; crosses - centroids in 30° cones around the p-axis; x-signs - centroids in 30° cones 
around the t-axis; dashed line is for the Cauchy rotation PDF with k = 0.075; left solid line 
is for the Maxwell rotation PDF dUD with cr$ = 7.5°. 
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Figure 11: 



Probability density function for rotation axes distribution. 
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GCMT catalog, N = 1 8, = 5.8, = 1 09.5 
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Figure 12: 

Distributions of rotation poles for pairs of focal mechanisms of shallow earthquakes in the 
GCMT catalog. Centroids are separated by distances between 0-50 km; magnitude threshold 
mt = 5.8; the rotation angle 109.5° < $ < 120°. 
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Simulation, Iseed = 2001 , N = 25000, NO = 1 754, zl = -0.05, zu = 0.05 
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Figure 13: 



Orientation distribution in the fundamental zone of the Rodrigues space for randomly rotated 
DCA sources. The points are shown in the central section of the fundamental zone 0.05 > 
X3 > —0.05 (see Fig. H]). The total point number is 25,000; 1,754 points are in the central 
section. 
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GCMT 1977-2012, = 5.8, N = 19397, NO = 8725, zl = -0.05, zu = 0.05 
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Figure 14: 

Orientation distribution in the fundamental zone of the Rodrigues space for shallow earth- 
quakes in the GCMT catalog. Centroids are separated by distances between 0-50 km; mag- 
nitude threshold = 5.8; the total number of events is 6,160. The points are shown in the 
central section of the fundamental zone 0.05 > X3 > —0.05 (see Fig. H]). The total number 
of earthquake pairs is 19,397; 8,725 pair points are in the central section. 
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